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Non-quasiparticle states in a half-metallic ferromagnet with antiferromagnetic s — d(f) 

interaction 


V. Yu. IrklihQ 

Institute of Metal Physics, 620990 Ekaterinburg, Russia 

Non-quasiparticle (incoherent) states which play an important role in the electronic structure 
of half-metallic ferromagnets (HMF) are investigated consistently in the case of antiferromagnetic 
s — d(f) exchange interaction. Their appropriate description in the limit of strong correlations 
requires a rearrangement of perturbation series in comparison with the usual Dyson equation. This 
consideration provides a solution of the Kondo problem in the HMF case and can be important for 
first-principle HMF calculations performed earlier for ferromagnetic s — d(f) interaction. 
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I. INTRODUCTION 

The s-d(f) exchange (Kondo lattice) model is widely used to describe many-electron (correlation) effects in d — and 
/—metals and their compounds. This model at large |/| is equivalent to the t — J model which is applicable, e.g., for 
high-Tc copper oxides. 

Besides the Kondo effect in the paramagnetic state, the correlation effects are important for magnetic semiconductors 
and half-metallic ferromagnets (HMF). The latter systems are widely investigated last time [1]. In particular, the 
interest in HFM is connected with possible applications in spintronics, including giant magnetoresistance (GMR) 
owing to spin-polarized electron transport across a non-magnetic metallic layer [2| (this phenomenon in HMF was 
predicted in Ref.[||). 

In the simple one-band case, HFM are strong (saturated) ferromagnets where (as well as in degenerate ferromagnetic 
semiconductors) the spin splitting exceeds the Fermi level, so that one spin subband is empty in the Hartree-Fock 
(Stoner) picture. However, this picture is radically changed by correlation effects which have in HFM a rather 
specific form, being connected with occurrence of incoherent non-quasiparticle (NQP) states flj. These states yield 
appreciable contributions to various physical properties (spin polarization, specific heat, transport characteristics etc.) 
pj, Hi. They occur due to quantum effects and are the only states with one of spin projections (cr =/ for the positive 
s — d(f) exchange parameter I and a =f for I < 0) near the Fermi level. The NQP states were also treated within 
first-principle calculations 

In such a situation, the s — d(f) exchange model provides an example of a non-trivial many-particle problem which 
can be investigated rather rigorously beyond perturbation theory. This is owing to that the ground ferromagnetic 
state is known exactly. As for field-theoretical methods, the case of a half-metallic ferromagnet can be investigated 
rather strictly. The exact solution in the one-electron case (empty conduction band) reduces the problem to the 
Schroedinger equation (two-body electron-magnon scattering problem). At finite temperatures (in the spin-wave 
region) and electron concentrations, a consistent consideration can be performed by constructing expansion in the 
occupation numbers of electrons and magnons. To this end, there can be used methods of equations of motion for the 
Green’s functions [§}, of the generating functional 0, and the expansion of the operator of evolution fo]. In this case, 
in each order there arise integral equations that describe the electron-magnon scattering (in their structure, they are 
similar to the integral equations of the Nagaoka type M)- 

However, the problem of NQP states in the case of negative (antiferromagnetic) s — d(f) parameter is still not 
solved satisfactorily, the physical picture in this case being especially non-trivial since the states below the Fermi 
become depolarized owing to occupied NQP states. In fact, the above-discussed expansion is somewhat different for 
ferromagnetic and antiferromagnetic s — d(f) exchange parameters (the difference is similar to that in the Kondo 
problem in the paramagnetic situation). 

In the present paper we present a consistent solution of the negative / problem and treat in detail this difference. 
In Sect.2 we derive the expressions for the one-electron Green’s functions, which are specific for the NQP situation. 
In Sect.3 we present the numerical calculation results and discuss them. Appendix deals with the limit of strong 
correlations. 
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II. CALCULATION OF ONE-ELECTRON GREEN’S FUNCTIONS 

We write down the Hamiltonian of the s — / exchange model in the form 

H — ^ -p ^ ^ '-fqS — qSq I ^^(SiCT a(7 / )o^ cr Cj cr ' , (1) 

k<r q icrcj' 

where cr are the Pauli matrices, fk is the bare band energy, I is the s — / exchange parameter. We consider the 
one-electron retarded Zubarev Green’s function. 


Gk a{E) — ((Cker I^kcr ))E 


( 2 ) 


We write down the Dyson equation for a ferromagnet at zero temperature 

G kt (£) = [E-t^-Z^iE)}- 1 , (3) 

where fkcr = tk — crIS is the mean-field spectrum, £k <r{E) is the self-energy. The latter quantity can be expressed as 
an irreducible Green’s function [§] 

E kt (£;) = (([Ck,Fint]| J ffint,c3 £ ]))| r (4) 

where Hi nt is the perturbation s — / Hamiltonian with the mean-field part being subtracted, “irr” means that the 
terms diverging as (E — t^ a )~ n should be omitted in the following expansion. 

To second order in / we have 




A± Oil 


7 k-ql. 


^ E — tk-qj. 

q 




Rt(E) = 


1 -riu- 


■k-qt 


^ E — tk- 

q 


■qt 


— u 


( 5 ) 

( 6 ) 


with ?i ko . = /(tk<r) the Fermi function, w q the magnon frequency (this is small as compared to characteristic electron 
energies, but provides a cutoff for logarithmic Kondo-like divergences). 

Performing decoupling at the next step (which corresponds to summing up the ladder diagrams 0 ) yields 


K(E) 


2J 2 Si% 
1 + o’/i? k 


( 7 ) 


For I > 0 we have in the half-metallic state (where = 0) 


G kt (£) = G^{E) = (E-tk + IS)- 1 

G^E) = G£ i(E ) = ( E - tk + ,S - 1 _;^ (£) ) ‘ (8) 

Thus, the spin-up electrons move freely and spin-down states are incoherent (the spectral density below the Fermi 
level comes from branch cut of the Green’s function owing to the resolvent R\^{E), but not from poles). 

The result © is exact in the case of empty conduction band since this corresponds to the exact solution of the 
two-body electron-magnon scattering problem. In the limit I —>• +oo Eq. © yields correct result (l23l) also for finite 
electron density. Thus for I > 0 the approximation m is quite satisfactory. 

In the limit I —>• —oo Eq.© gives correctly the spectrum of spin-down quasiparticles (see Ea. (l23l) h However, for 
I < 0 the ladder approximation is not appropriate for describing NQP states, although it reproduces a non-pole 
structure. There are several problems. First, Eq. © does not reproduce atomic limit and limit of strong correlations 
(see Appendix). Second, it contains the false Kondo pole, which may lead to some problems. 

Thus we use, instead of ©, another representation for the Green’s function, which is appropriate for describing 
NQP states. Writing down the equations of motion we obtain 

G kt (£) = G°^(E) + I 2 [Gi t (E)]H k (E) 

$k (E) = 5> kp (E) 


( 9 ) 
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with 

r k p (E) = Y / ((SpC k - pl + ss;c k . pt \ci +ri s+ + cl +rt ss z r )) E 

pr 

= X]^ Sf P C k-p4-l C k+r4. 5 r + 2cj c + r t <5S r»B ( 10 ) 

pr 

Here 5A = A — (A), we have used symmetry, purely longitudinal fluctuation terms (without spin flips) being not 
important at low temperatures. Note that the Green’s function (flTill is not irreducible. The first term in ([9]) does not 
give a contribution to the density of states in the HMF state and, in any case, for large |/|. 

The equation of motion for T reads after decoupling 

(E — fk-p — IS) T kp (E) = 2S'rik-p4, — 2 ^~^(c k _ at c k - p j - 5 , p _ c| ) 

q 

-2J5n k _p^((c kt |ct +r; 5+ + 24 +Pt 5S r *)) B (11) 


where n kcr = (c ko .c kcr ). Using the equation 

(E - t k + IS)J2Mcl +ri S+ + 2ct +rf 6S z )) E = -/£ T kp (£) (12) 

r p 

we derive the integral equation 

(E — f k - P — IS) r kp (S) = 2Sn,k-pi(E) — I^ ^ ] T kq (l?) (13) 

q 

with the renormalized occupation numbers 

hki = n n - - 53^k-qt c H‘ S 'q') ( 14 ) 

q 

The second term in Q is somewhat similar to that in the paramagnetic Kondo problem (e.g., in the Nagaoka 
decoupling 0 ). Calculating the corresponding Green’s function 

((S+Ck-qtlc^))^ = -2 ISG^(E)G 0 n (E) (15) 

and using the spectral representation we obtain 

( C k -qt Ck 4-'-’q ) = ~ IlSn^l/ (A k | — t k _ q + IS) (16) 

The spin-down distribution function can be calculated as 


where the residue is given by 


«14 = J(E^) 


Z kl =(l-—Xi(E = E n ) 


(17) 


At large |/| we can use the simple approximation (see Appendix) 


7 25 -v ( 2S + 1 \ _25+l , 

Zki ‘2S -4- 1’ ^^1 \ oc I q J \Ek\.) 


2 S 


2 S 


(18) 


It should be noted that these calculations are not quite strict, but are in fact in spirit of the quasiclassical 1/S 
expansion Q, treatment of higher orders being needed for justification. 
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FIG. 1: Non-quasiparticle spin-up density of states N(E) in a half-metallic ferromagnet for the semielliptic bare band with the 
width of IF = 2 and different I values, S' = 1/2 (a) and S = 3/2 (b) according to (1211) 


Ea. lfldl) can be easily solved. In higher orders of perturbation theory (taking into account next terms in the 
expansion) the denominators (E — l k _ p — IS)~ are replaced by the exact Green’s functions G£ j_(.E), so that we 
derive 

d> (E) = „ 

1 + IRl(E)G° kt (E)/Gi l (E) 

Rl(E) = £n k _^(2?)G^(i5 + u> q ) 

q 

This structure differs from the “ladder” one 0 by the presence of the factor G^(E)/G^(E). 


(19) 

( 20 ) 


III. RESULTS OF NUMERICAL CALCULATIONS 


The final result for the NQP contribution to the spin-up Green’s function can be written down in the form 


G kt (E)^G° t (E) 


_2 ISGj t (E) _ 

[G^GE)]- 1 + [IRl(E)Gl t (E)]-i 


( 21 ) 


Unlike 0, expression m has no false “Kondo” pole, so that analytical properties are fair. Besides that, it yields 
the correct structure in the large-|/| limit (see Appendix). 

In our model, HMF state is characterized by large spin splitting which exceeds the Fermi energy. Thus the usual 
band states are absent below the Fermi level, and the NQP states owing to branch cut of the resolvent (1^01) dominate 
in this region. The instability of HMF state is determined by occurrence of the pole of the Green’s function m 
below Ep, which can take place with increasing the Fermi energy. 

Despite the presence of the energy gap in the spin-up subband, the current carriers are depolarized even in the 
ground state owing to the electron-nragnon interaction. Physically, the occupied spin-down electron states are a 
superposition of the states |<S)|4.) and \S — 1)| t)> so that we can pick up from these states spin-down and spin up 
electrons with the corresponding spectral weights. The depolarization is maximum in the large-|/| limit where (for 
S = 1/2) it reaches 100%. 

Fig.l shows the energy dependence of NQP density of states for different values of the s — f exchange parameter, 
the occupation numbers being taken in the approximation (1181) . The energy is calculated from the centre of the bare 
band shifted by atomic value I(S + 1). The chemical potential corresponds to — 1 < t k — I(S + 1) < /To = —0.7. Note 
the shift of the true chemical potential /r with the change of I owing to spin-polaron renormalization. One can also 
see that the |/| dependence of the NQP contribution size is weak since \I\ is in fact not a small parameter in the HMF 
regime (large spin splitting). 
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FIG. 2: Non-quasiparticle spin-up density of states in a half-metallic ferromagnet with different I values S = 1/2 (a) and 
S = 3/2 (b) according to 0 


For comparison, the results according to Eq. © are presented in Fig.2. The chemical potential corresponds to 
— 1 < tk ^ IS < Ho = —0.7. One can see that the visual picture is quite different, especially in the quantum (small 
5) case. In particular, this approximation does not yield band narrowing (which is 25/(25' + 1) in the large-/ case). 
A complicated behavior at intermediate I for 5 = 1/2 demonstrates the influence of the Kondo singularity. 

Note that in fact the Kondo divergence has no crucial effect because of the one-side character of the singularity, so 
that N^(Ep) = 0; besides that, smearing due to spin dynamics should be taken into account. Thus a dependence on 
cutoff at the spin-wave frequency arises, which is stronger in the approximation ©. 


IV. DISCUSSION AND CONCLUSIONS 

The problem of half-metallic ferromagnetism in the s — f exchange model is somewhat similar to the Kondo lattice 
problem in paramagnetic (or weakly magnetic) state, but can be investigated rigorously by rather simple method of 
the Green’s functions. 

In the present work, we have obtained a consistent solution of the Kondo problem in the case of strong (half-metallic) 
ferromagnetism in terms of NQP states, the ultimately strong correlation limit being correctly described. In our case, 
there are no Kondo corrections to the ground state magnetic moment, so that ferromagnetism remains saturated. 
On the other hand, in the case of weak magnetism the standard scaling theory starting from the approximation © 
describes successfully the competition of the Kondo effect (screening of magnetic moments) and intersite exchange 
interactions El- Note, however, that half-metallic solutions with a Kondo-reduced magnetic moment can occur in 
the strong-coupling regime [Hi. 

The Kondo divergences in a ferromagnetic state were also considered for the usual magnetic [j~4j] and pseudospin 
El problems, transport properties (e.g., sharp energy dependence of resistivity and large thermoelectric power in the 
case of impurity scattering) being discussed. 

The HMF spectrum picture obtained includes to first line non-quasiparticle states which are the only electron states 
near the Fermi level for one of spin projections. The example of HMF turns out here instructive to investigate a very 
important problem of passing from one-electron to many-electron statistics. 

The NQP states were classified as minority ones (above the energy gap, e.g., in NiMnSb and CrC >2 1) and majority 
ones (below the energy gap, e.g., in Mn 2 VAl @). The latter situation just corresponds to the case I < 0. Thus our 
treatment can be used in a combined consideration with first-principle band calculations. 

The author is grateful to M. I. Katsnelson and A. I. Lichtenstein for useful discussions. This work was supported 
in part by the Division of Physical Sciences and Ural Branch of Russian Academy of Sciences (project no. 15-8-2-9). 
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Appendix. Atomic and strong correlation limits 


a 


The atomic-limit Green’s functions can be calculated exactly for arbitrary electron density n and magnetization 
. They have a two-pole structure corresponding to two atomic levels 


G G (E) 


1 / S + n + cr(S) 

25+1 V E +IS 


S + 1 — <r(5) — 7l\ 
E-I(S+ 1) ) 


( 22 ) 


Note that a part of spectral weight vanishes in the large-1/| limit. 

The strong correlation limit \I\ —> oo can be also investigated rigorously [16}. The solution can be found by using 
the atomic representation of Hubbard’s X-operators. (However, for finite |/| such analytical calculations in a closed 
form are difficult and not instructive: there are too many states with different spin projections and occupations.) 

In the limit / —> +oo we have 


Gl(E) 

^k(e) 


e - tk 


. G l k (E) = 


e-tk 


2 S 

R^)\ 


-1 


1 fit k—q) 

q 6 " ^ ~ ^ 


with e = E + IS. 
In the limit I — 


—oo we derive 


o q 

Gtw = 2 STT^)' 


(23) 

(24) 


(25) 


Gl(E) = 


2 S 


2S + 1 


e — tt + 


2 S 


KM! 


K(e) 


V k-q) 

^ e-tt 


‘'k-q 


(26) 

(27) 


with e = E — I{S + 1), t k = [25/(25 + l)]tk- 

The Green’s function (l26l) has no poles, at least for small current carrier concentration, and the whole spectral 
weight of minority states is provided by the branch cut (non-quasiparticle states). The one-electron distribution 
functions read 


^kj. 


25 

25+1 


f(t k), n kt 


n 

25+1’ 


(28) 


so that rikf is very weakly k-dependent. One can see that the equation for the chemical potential, ^ k(7 n k<r = n, is 
satisfied just owing to NQP spin-up states. 

It should be noted that the motion of observable strongly correlated quasiparticles are described by many-electron 
projection X-operators, 


V25TT^’ 4 = M-l/2,5-l/2)(i,5,5| 

where \i,S',M) is the state on a site i with total spin S' and its projection M. The corresponding distribution 
function does not contain the factor of 25/(25 + 1). 


((a+l+L))^ = k)> (4 x k) = fit k)- 


( 29 ) 
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